Thermodynamics of the glassy state: effective temperature 
as an additional system parameter 
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A system is glassy when the observation time is much smaller than the equilibration time. 
A unifying thermodynamic picture of the glassy state is presented. Slow configurational modes 
are in quasi-equilibrium at an effective temperature. It enters thermodynamic relations with the 
configurational entropy as conjugate variable. Slow fluctuations contribute to susceptibilities via 
quasi-equilibrium relations, while there is also a configurational term. Fluctuation-dissipation rela- 
tions also involve the effective temperature. Fluctuations in the energy are non-universal, however. 
The picture is supported by analytically solving the dynamics of a toy model. 
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Thermodynamics is an old but powerful subject. Af- 
ter the invention of the steam machine it was needed 
for understanding their optimal efficiency. It applies to a 
wide variety of systems, ranging from ideal gases to black 
holes. Equilibrium thermodynamics is a well understood 
subject. Important contributions were made by Carnot, 
Clausius, Kelvin, Boltzmann, and Gibbs. 

Thermodynamics for systems close to equilibrium was 
worked out in the mid of this century. Applications are 
systems with heat flows, electric currents, and chemical 
reactions. A basic assumption is the existence of local 
thermodynamic equilibrium, and the task is to calculate 
the entropy production. Important contributions were 
made by de Donder, Prigogine, de Groot and Mazur. 

Thermodynamics for systems far from equilibrium has 
long been a field of contradiction and confusion. A typ- 
ical application is a window glass. A cubic micron of 
glass is not a crystal, it is an undercooled liquid which 
has fallen out of equilibrium. In a glassy system the re- 
laxation time T eq (T) of the slow or so-called a-processes 
has become larger than the observation time, while the 
fast or /3-processes are still in equilibrium. Waiting long 
enough might bring the system back to its equilibrium 
state. In glasses one often assumes the Vogel-Tamman- 
Fulcher law r eq ~ exp(A/ (T — Tk))- In many other sys- 
tems one encounters an Arrhenius law r eq ~ exp(^4/T). 
As we shall demonstrate on a toy model, such systems 
display the same glassy behavior. 

We shall consider glassy transitions for liquids as well 
as for random magnets. The results map onto each other 
by interchanging volume V, pressure p, compressibility 
K = —d In V/ dp, and expansivity a — d In V/ dT , by mag- 
netization M, field H, susceptibility x — (l/N)dM/dH, 
and "magnetizability" a = (—l/N)dM/dT, respectively. 

A glass forming liquid exhibits near the glass tran- 
sition smeared discontinuities in quantities such as the 
heat capacity C p , the expansivity and the compressibil- 
ity. This defines a smeared glass transition line T g (jp) or 
p g (T), with behavior similar to continuous phase transi- 



tions of classical type, i.e., with specific heat exponent 
a = 0. Since the twenties quite some attention has 
been payed to introduce a thermodynamic description, 
see e.g. || ||. It was investigated, in particular, 
whether the discontinuities satisfy the two Ehrenfcst re- 
lations (the analogs for second order transitions of the 
Clausius-Clapeyron relation of a first order transition), 
and whether the Prigogine-Defay ratio (see eq. (|25|)) 
equals unity. Very recently we have explained the ex- 
perimental observations 0. 

A state that slowly relaxes to equilibrium is charac- 
terized by the time elapsed so far, the "age" or "waiting 
time" . For glassy systems this is of special relevance. In 
spin glass experiments non-trivial cooling-heating cycles 
can be described by an effective age || . One thus char- 
acterizes a non-equilibrium state by three parameters, T, 
H or p, and the age t or the cooling rate T = dT/dt. For 
thermodynamics a more suitable third variable is the ef- 
fective temperature T e {t), introduced half a century ago 
by Tool §} @. T e ^ T will describe the best quasi- 
equilibrium the a-processes could reach so far; it follows 
in simple cases by equating T eq (T e ) = t. 

For a set of smoothly related cooling experiments T, it) 
at fields Hi one may express the effective temperature as 
a continuous function: T e/ i(t) — > T e (T,H). This sets a 
surface in (T,T e ,H) space, that becomes multi-valued if 
one first cools, and then heats. To map out the whole 
space many sets of experiments are needed, e.g., at dif- 
ferent cooling rates. Thermodynamics amounts to give, 
for a certain class of systems, universal relations between 
state variables at nearby points in this space. 

Within this framework some new results were obtained 
for a spin glass model with one step of replica symme- 
try breaking || || such models are related to sys- 
tems without disorder that have one, and only one, di- 
verging time scale. The fast and slow modes do not 
only have their own temperature, they also have their 
own entropy: the entropy of "equilibrium processes" 
S cp , and the "configurational" or "information" entropy 
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or "complexity" X , respectively. The total entropy is 
S = S cp + T. Previous results can be summarized in 
setting dQ = TdS cp + T e dl. In the coarse of time the 
system will satisfy (T — T e )dl > 0, as required by the 
second law. In combination with the first law, the ther- 
modynamic relations can then be represented as 



dU = TdS cp + T e dl +dW 
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U(T(t),t,H) -> U(T,T e {T,H),H) and M is, of course, 
not satisfied; however, the above implies that its violation 
is related to the change of configurational entropy: 
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Here we wish to add to this scheme the non-equilibrium 
susceptibilities. We will derive the results in a toy model 
with fields H a and magnetizations M a (a — 1,2). As 
could be guessed from M a = M a (T(t),T e (t, H), H) with 
T e (t, H) — * T e (T, H), there appear two terms: 
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First there is the expected fluctuation contribution 



changes much faster than T e (t'), allowing to replace T e (t') 
by T e (t) in the aging regime, and by T in the equilibrium 
regime t' « t. 

Note that the ratio dt'C/G depends on time. The sit- 
uation with constant T e is well known from mean field 
spin glasses jl3). However, this time-independence is an 
artifact of the mean field approximation || . 

Finally, one estimates scaling in the aging regime of 
two-time quantities as 



C(M')«C(- 
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showing immediately the familiar t/t' scaling; there may 
be logarithmic scaling corrections. 

The full picture (eqs. (l)-(9)) says that slow modes are 
at a quasi-equilibrium at T e . As in a plasma, slow and 
fast modes equilibrate at their own temperature. 

In a number of simple models there occurs a dynami- 
cal glassy state when cooling near T — 0. If there is only 
one time-scale, there remain only slow processes in the 
frozen phase (only a-, no /3-processes). This implies that 
U has no explicit dependence on T, and that S op = 0. 

In order to corroborate our statements, we consider a 
toy model involving free spherical spins (%2Sf = N). 
The Hamiltonian contains two parts, a "self- interaction" 
term involving quenched random fields Tj = ±r with 
average zero, and a coupling to an external field H 
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Notice that fast and slow processes again enter with their 
own temperature, while the third term is small. Since 
T e ^ T, there occurs a new, configurational term 
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It originates from the difference in the system's structure 
for cooling experiments at nearby fields Jll| . Such uni- 
versal relations would not hold for energy fluctuations. 

We shall also see in the toy model, that the correla- 
tion function C ab (t, t')=(l/N)(5M a (t)SM b (t')) and the 
response function G a b(t,t')=(l/N)SM a (t)/SH b (t') sat- 
isfy the fluctuation-dissipation relation 



dC ab (t,t') 

dt> 



T e (t')G ab (t,t') 
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In terms of the "staggered" magnetization M s = 
(1/r) J2i T *Si on e simply has H = -TM S - HM. 

In equilibrium at low T the internal energy reads 
U eq /N = -K + ±T, the entropy S eq /N = \ \xi(eT/K), 
the magnetizations M eq /N = H/K - \HT/K 2 , and 
M s , eq /N = T/K - ±TT/K 2 , where K = VT 2 + H 2 . 



This model gets glassy behavior when it is subject to 
Monte Carlo dynamics with parallel update, which cou- 
ples the spins dynamically. The time evolution can be 
solved exactly, since it maps closely on the dynamics for 
uncoupled harmonic oscillators, introduced by Bonilla, 
Padilla, and Ritort 14 1, when extended to include a field. 

Per time step 1/N one makes parallel Monte Carlo 
moves, Si — ► S[ = Si + ri/yN, with Gaussian noise 
having <ri>~ and <r 2 >= A 2 . Next one makes a 
global rescaling of the length of the spins in order to 
keep S^ 2 — N. This leads to the final update 



Sa — Si 
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in the aging regime, while in the equilibrium or short- 
time regime T, replaces T e . This been confirmed numer- 
ically for a soft sphere glass Eq. (^) is consistent 
with (||): in performing L G ab (t, t')dt' one uses that dfC 



Denoting TL by E, it is simple to calculate the joint tran- 
sition probability of x = E' — E and y = M' — M: 



p(x,y\E,M)=p(x\E)p(y\x,E,M) 



(12) 
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where both factors are Gaussian. In terms of 

M H H 



N ' 



K + KS* (13) 



the centers of the Gaussians are 



x = \A 2 {K-e): y = 
while their variances read 
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respectively. Following Metropolis, a move E' = E + x, 
M' = M + y is accepted with probability W((3x), where 
W = 1 if x < 0, while W = exp(-/3x) if x > 0. The dy- 
namics is now fully specified. The average energy evolves 
according to a closed equation for s(t) Q, 



de 
dt 



dxW (f3x)xp(x\e) 
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We find that the average magnetization satisfies similarly 

— = / dxW((3x)y oP (x\s) (17) 
at J^oo 

At low T there occur Arrhcnius laws for the equilib- 
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rium relaxation times viz. 
{irT/MA) 1 / 2 exp(A/T) with A = A 2 K/8, responsible for 
non-equilibrium behavior in typical cooling procedures. 

We first consider the situation at T = 0. Starting from 
a random initial condition, the system will slowly evolve 
towards the ground state. It is this evolution that we 
now wish to capture within a thermodynamic framework. 
Equating U(t) = U eq (T e ,H) leads to Q 



Te(t) 



(18) 



As asserted above, the same T e is obtained via the relax- 
ation time: * = T eq {f e {t)) -> f e (t) = T e (t) + 0(T 2 ). 

In equilibrium it holds that in = 0. The result m ~ 
1/t r*> exp(— (3 e A) proves that the magnetization very 
closely follows its quasi-equilibrium value m eq {T e ,H). 

To test the thermodynamics we need the configura- 
tional entropy. It is defined as the logarithm of the num- 
ber of states leading to U(t,H) — U eq (T e , H). Since 
Sep = it simply holds that X = S eq (T e (t), H). We 
can now verify the relation dU — T e dl — MdH ~ Af s dr. 
At constant H and V it is valid, because one has replaced 
T — > T e in energy and entropy. Then one can take the dif- 
ference of two evolution experiments at two nearby H's. 
The relation remains satisfied since M(t) sa M eq (T e , H). 
The modified Maxwell relation (Q) is also obeyed. Fi- 
nally one can change H in the coarse of time. Then 
m < HT e /K 2 as long dH(t)/dT e {t) < [3 e (t)H{t)A 2 . 
This is a mild condition. If H is changed quicker, or if it 



goes to zero too rapidly, the system will not be able to 
reach a quasi-equilibrium described by T e alone. 

We have also considered the fluctuation formula. They 
are too lengthy to be reproduced here. We found 
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The latter two relations imply that T 2 {5M 2 ) w H 2 (SM 2 } 
w — TH(SMSM S ). Since there are no fast processes, it 
immediately follows that the quasi-equilibrium relation 
(g) is satisfied for all four cases, to leading order in T e . 
In contrast, for M a — > Tt, H a — > j3 or j3 e the correspond- 
ing relations are violated. 

The two-time correlation functions C a b(t,t') and re- 
sponse functions G a b{t,t') exhibit related behavior. 
These functions decay as C a i,(t,t') = C a b(t' ,t')h(t')/h(t) 
with h(t) = exp(f3 e (t)A) ~ Whit and satisfy the 
fluctuation-dissipation relation (^) for any non-zero value 
of H and Y ijTsf . (For the oscillator model at zero field 
these relations were derived in p4|.) 

Now we consider a cooling experiment from high tem- 
peratures. A glass transition will occur when the cooling 
timescale T/\T\ becomes comparable to the equilibration 
timescale T eq (T e (t)). Let us assume that it happens at a 
low temperature T g <C A, so at an exponential time scale 
t g = T eq (T g ) ~ expf3 g A. This will imply that the width 
of the transition region is small: AT g ~ T 2 /A <C T g . As- 
suming to have in that region a constant dimensionless 
cooling rate Q = TT' eq (T), we derive from ( |l6| ) 
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This equation is similar to but different from earlier pro- 
posals H [ [Tcfl , and might be universal for narrow glassy 
transitions. It covers three cases: a) Q > 1: normal cool- 
ing towards or in the glass; b) < Q < 1: cooling in a 
glassy state so slowly that equilibrium is achieved later; 
c) Q < 0: heating up in the glassy state. The apparent 
specific heat c = ^dT e /dT decreases if cooling is not too 
slow. In a heating experiment it is negative, but it pro- 
duces near T g the well know overshoot c > h , with height 
and shape depending solely on Q. 

In a cooling experiment from large T eq. (|20|) yields 
initially T e (t) = T(t), up to exponential corrections, de- 
scribing thermodynamic equilibrium at the instantaneous 
temperature. Below the glass transition region one has 
T e nT g + (T -T g )/Q « T- q x (t). This agrees with @ 
and shows that the actual temperature and the cooling- 
history are irrelevant: to leading order the energy just 
evolves as if the system had been quenched to zero tem- 
perature, and aged there. It can be checked that all re- 
lations linear in T e remain as at T = 0, thus supporting 
the picture (l)-(9) . 
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This solution allows us to check the Ehrenfest rela- 
tions along the glass transition line H g (T) in cooling pro- 
cedures (Q > 1). Below the transition region one has 
c = 1/(2Q), a = H/(2K 2 Q) and X = T 2 /K 3 - (H/2K 2 ) 
dT e /dH. Comparing with the paramagnet and using 
that along the transition line the equality T e = T implies 
that dT e /dT+(dT e /dH)(dH g /dT) = 1, we find that the 
jumps in a and \ satisfy the first Ehrenfest relation 



Aa = A X -^ 



(21) 



Until very recently, it was widely believed that this rela- 
tion is violated at the glass transition. We pointed out 
that experimentalists had inserted some short-time value 
of x or k II, like the "zero-field-cooled" susceptibility 
Xzfc, he. the first term of ^ fluct in eq. (||). In spin 
glasses with one step replica symmetry breaking xzfc 
is strictly smaller than the X °f the paramagnet [fl6|| . 
In the present model, there are no /3-processes, so Xzfc 
even vanishes. On the other hand, the long-time or field- 
cooled value XFC=X fiuct is continuous at T g . The cor- 
rect discontinuity, A\ = — x conf arises from eq. (0). As 
Aa = (dm / dT e )(dT e / dT — 1), this explains in detail why 
the first Ehrenfest relation is satisfied automatically Q . 

It was also shown in Q that the modified Maxwell re- 
lation (||) leads to the modified second Ehrenfest relation 
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We can now verify that it is also satisfied. The new last 
term needed for validity beyond equilibrium. 
The Prigogine-Defay ratio may be expressed as 

ACA X 1 dT dl 

~ NTg(Aa) 2 ~ NAa dT > dH { 6) 

The definition looks as an equilibrium relation, and it 
shown that II must be larger than unity for mechani- 
cal stability 0. This was based, however, on the in- 
valid assumption of thermodynamic freezing of a set of 
unspecified order parameters. The equivalent relation 
II = (AC / NT g Aa)dT g / dH allows II < 1 if dT g /dH, 
which depends on dQ/dH, is small enough. We real- 
ized that already in the classic experiment of Rehage and 
Oels |l7j on the glass transition in atactic polystyrene 
there occurs a value II = 0.77 In our present model 
IT = (K 2 / HTg)dT g / dH becomes less than unity when- 
ever dXg/dH is positive. This condition occurs in half of 
the smoothly related sets of cooling trajectories. 

In conclusion, we have proposed a unifying thermody- 
namic picture of the glassy state. It does not apply to 
ideally slow experiments, but to conditions that are typ- 
ically met. As in a plasma, slow and fast modes equili- 
brate at their own temperature. On long time scales and 
under mild conditions, global thermodynamical quanti- 
ties, like the energy and volume or magnetization(s), go 
to quasi-equilibrium values at a certain effective tempera- 
ture T e (t). Slow fluctuations contribute to susceptibilities 



with factor 1/T e . For cooling trajectories at two nearby 
external fields (pressures), the difference of the magne- 
tization (volume) involves the usual fluctuation suscep- 
tibility (compressibility), and a new structural contribu- 
tion, that arises from the difference in effective temper- 
atures. Correlation and response functions exhibit t/t' 
scaling, and satisfy a fluctuation-dissipation relation in- 
volving T e . 

We have verified our picture in the glassy phase of 
the low-temperature dynamics of an exactly solvable toy 
model, that contains two external fields. We expect that 
the fields may also stand for other "mechanical" forces, 
such as pressure or chemical potential. 

The picture has the right signs to be valid for a class 
of glassy systems. For verifying it in glass forming liq- 
uids it is desirable to map out the full (T, T e ,p)-space 
by cooling experiments, and to check that against aging 
experiments. 
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